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Land surface vegetation phenology is an efficient bio-indicator for monitoring ecosys- 
tem variation in response to changes in climatic factors. The primary objective of the 
current article is to examine the utility of the daily MODIS 500 m reflectance anisotropy 
direct broadcast (DB) product for monitoring the evolution of vegetation phenological 
trends over selected crop, orchard, and forest regions. Although numerous model-fitted 
satellite data have been widely used to assess the spatio-temporal distribution of land 
surface phenological patterns to understand phenological process and phenomena, cur- 
rent efforts to investigate the details of phenological trends, especially for natural 
phenological variations that occur on short time scales, are less well served by remote 
sensing challenges and lack of anisotropy correction in satellite data sources. The daily 
MODIS 500 m reflectance anisotropy product is employed to retrieve daily vegetation 
indices (VI) of a 1 year period for an almond orchard in California and for a winter wheat 
field in northeast China, as well as a 2 year period for a deciduous forest region in New 
Hampshire, USA. Compared with the ground records from these regions, the VI trajec- 
tories derived from the cloud-free and atmospherically corrected MODIS Nadir BRDF 
(bidirectional reflectance distribution function) adjusted reflectance (NBAR) capture 
not only the detailed footprint and principal attributes of the phenological events (such 
as flowering and blooming) but also the substantial inter-annual variability. This study 
demonstrates the utility of the daily 500 m MODIS reflectance anisotropy DB prod- 
uct to provide daily VI for monitoring and detecting changes of the natural vegetation 
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phenology as exemplified by study regions comprising winter wheat, almond trees, and 
deciduous forest. 


1. Introduction 

Land surface phenology, defined as the seasonal pattern of variation in vegetated land 
surfaces characterized by satellite remote sensing products (Jolly, Nemani, and Running 
2005; Zhang, Friedl, and Schaaf 2006), is an efficient bio-indicator for monitoring terres- 
trial ecosystem variation in response to changing climatic factors. Changes in phenological 
events are increasingly being seen as important signals of year-to-year climatic variations 
or even global environmental changes, and phenological studies are being used to moni- 
tor the impact of climate change on individual species at a number of scales (Reed et al. 
1994; Schwartz, Ahas, and Aasa 2006; Cleland et al. 2007). The timing of leaf development 
in many temperate deciduous species correlates with cumulative springtime temperatures 
(Canned and Smith 1983; Rotzer, Grote, and Pretzsch 2004; Zhang and Goldberg 2010; 
Zhang, Goldberg, and Yu 2012). The length of growing season may strongly affect the 
associated seasonality of ecosystem-scale carbon exchange, especially in deciduous forests 
(White, Running, and Thornton 1999; White et al. 2009; Richardson, Braswell, et al. 2009, 
2010, 2012). A number of studies have linked vegetation’s response via phenology and 
length of regional growing seasons with climatic and ecological changes (e.g. Myneni et al. 
1997; Parmesan and Yohe 2003; Liang and Schwartz 2009; Liang, Schwartz, and Fei 2012). 
Conversely, various studies show vegetation’s feedback on local climate via biochemical 
and biophysical processes. Temperature is the main factor controlling the seasonality of 
vegetation growth in humid temperate climates (Schwartz 1998; Zhang et al. 2004), and 
water availability governs vegetation phenology in arid and semiarid ecosystems and in 
seasonally dry tropical climates (Justiniano and Fredericksen 2000; Kramer, Leinomen, and 
Loustau 2000; Zhang et al. 2010) where rainfall activates the emergence of green leaves and 
controls the duration of vegetation growth (Spano et al. 1999; Cook and Heerdegen 2001). 
The pattern of precipitation during rainy seasons is critical to crop germination, growth, 
and harvest (e.g. Lobell and Field 2007; Osborne et al. 2007; Omotosho 1992), and the 
early or late arrival of precipitation events can affect estimations of agriculture growth and 
yield. 

Many land surface phenology investigations have been devoted to the development of 
various methods (de Beurs and Henebry 2010) for identifying and calculating phenolog- 
ical features using time series of vegetation indices (Vis). A regression-based model of 
first leaf emergence is used to detect the onset of spring (Schwartz 1990). User-defined 
thresholds depending on land cover types are adopted to identify the phenological dates 
(White, Thornton, and Running 1997; White et al. 2002; Schwartz, Reed, and White 2002). 
In most recent studies, phenological time series are explored using the TIMESAT smooth- 
ing algorithm (Jonsson and Eklundh 2004; Tan et al. 2011) to fill gaps and reduce noise 
in the satellite signal, through a complex six-parameter model to fit sparse greenness pro- 
files (Badhwar 1984), and also through piecewise logistic functions to determine onset 
dates (Zhang et al. 2003; Ganguly et al. 2010). These methods have significantly improved 
our understanding of spatial phenological patterns (Zhang and Goldberg 2010), the link 
between climate change and ecosystem variations (Busetto et al. 2010), and the contri- 
bution of forest phenology to the estimation of carbon flux (Jenkins et al. 2007). On the 
other hand, fitted or smoothed phenological time series over multi-day composite prod- 
ucts can mask high-frequency vegetation changes, possibly missing short-term changes 
evident in crop defoliation, post snowmelt green-up in the Arctic, or fire and other sudden 
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disturbances (White, Thornton, and Running 1997; White and Nemani 2006; McKellip 
et al. 2005; Narasimhan and Stow 2010; Ju et al. 2010). As some resource management 
actions rely on knowing the date of such disturbances as accurately as possible and react- 
ing to the disturbance as quickly as possible afterwards (e.g. replanting when crops fail or 
reducing erosion after fire), fine-temporal resolution satellite products are required for the 
exploration and monitoring of land surface vegetation phenological events. 

Traditionally, land surface phenological metrics have been extracted from weekly, 
10 day, and 14 day NOAA/AVHRR (Advanced Very High Resolution Radiometer) archives 
(Reed et al. 1994; Schwartz, Reed, and White 2002; de Beurs and Henebry 2008; Kross 
et al. 2011), and for the past decade from 16 day MODIS operational products (Zhang, 
Friedl, and Schaaf 2006; Ganguly et al. 2010). Compared with AVHRR, MODIS has 
seven carefully selected bands for land applications, finer spatial resolution (250 m to 
500 m), improved geolocation accuracy (50 m at nadir), systematic atmospheric correc- 
tion schemes and cloud screening, and upgraded sensor calibration (Justice et al. 1998; 
Wofle et al. 2002; Box et al. 2006; Heidinger, Cao, and Sullivan 2002). The routine 
MODIS reflectance anisotropy and albedo products (MCD43) (Schaaf et al. 2002, 2011) 
use high-quality, multi-spectral, cloud-free, atmospherically corrected surface reflectances 
from both Terra and Aqua to provide global quantities with the anisotropy correction for 
directional reflectance at a 500 m spatial resolution every 8 days (based on a 16 day win- 
dow). While the operational product was originally envisaged as a daily product, archive 
constraints restricted it to an 8 day retrieval schedule. However, due to the increasing needs 
of the MODIS direct broadcast (DB) user community, a daily DB algorithm was imple- 
mented (Shuai 2010) and released through the University of Wisconsin-Madison Space 
Science and Engineering Center (SSEC) to more than 125 MODIS DB stations around 
the world (http://cimss.ssec.wisc.edu/imapp/db_brdf_vLO.shtml). This package, provid- 
ing daily 500 m albedo measures and Nadir BRDF (bidirectional reflectance distribution 
function)-adjusted reflectance (NBAR), is well suited for investigating the natural trajec- 
tories of various vegetation species with rapid growth patterns or impacted by frequent 
disturbances. 

In this study, daily MODIS 500 m NBAR-based phenological markers are compared 
to ground measurements carried out over three homogenous sites (an almond orchard 
in California, a winter wheat region in northeast China, and a deciduous forest in New 
Hampshire, USA). More precisely, the primary objective is to explore the utility of daily 
MODIS-derived NBAR to characterize the VI trajectories of various vegetation species. 
A secondary objective is to examine the natural phenological changes in the above three 
vegetation species. A further objective is to maximize the inherent ability of the MODIS 
DB system by providing the near real-time phenological information on vegetated areas for 
the pressing requirements in rapid-response systems. 


2. Data and methods 
2.1. Field data 

The first study area (37° 30' 41.04" N, 121° 12' 47.73" W), located in western Stanislaus 
County in the northern San Joaquin valley region of central California, USA, is a typical 
almond orchard with more than 100 acres of continuous rows of mature trees. In such 
orchards, grass, or planted crops are often found in the understorey vegetation. Since 
2001, weather conditions and phenological events related to California almond orchards 
have been reported based on field observations from numerous growers through the 
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blue diamond growers website (http://www.bluediamond.com/applications/in-the-field). 
The regularly scheduled field reports describe the detailed status and evolution of each 
growth stage from dormant, green tips, buds, blooming, jacketing, to the nuts developing 
through descriptions and field photos. In this case study, we extracted the phenological 
dates (Table 1 ) in terms of the estimated percentage of individuals entering each stage for 
the 2006 crop year of the selected orchard. 

The ground-based crop growth stages for the second site were collected over the 
Yucheng Experiment Site (YCES) in the agricultural region of the northeast China plain by 
the Chinese ecosystem research network at YCES. This site, located southwest of Yucheng 
city at 36.8290°N, 1 16.5704°E, has been planted with winter wheat and summer corn over 
90% of its area since the year 2000. The in situ data record used in this case study includes 
frequent field specification of the surface phenology information during the 2005 growing 
season. Table 1 summarizes the ground observations of the phenological period available 
for the winter wheat crops at Yucheng. According to the field criterion, the ‘start’ of a 
growth stage means that ~ 10% of the individual crops have entered into the growth phase, 
while the ‘end’ indicates that "-90% crops have completed this stage. 

The third case study was conducted at the Bartlett Experimental Forest (BEF, 44° 17' 
N, 71° 3' W), located within the Saco Ranger District of the White Mountains National 
Forest in north-central New Hampshire, USA. The dominant trees are northern hardwood 
species ( Acer sacchamm, sugar maple; Fagus grandifolia, American beech; and Betula 
alleghaniensis, yellow birch) with some red maple, paper birch, and occasionally conifers, 
such as eastern hemlock (Richardson et al. 2007; Richardson, Flollinger, et al. 2009). 
A digital webcam mounted near the top of a 26.5 m tower, looking north, and angled 
slightly downward, has been providing digital images recorded daily between 12:00 and 
14:00 hours since 2005 (Jenkins and Richardson et al. 2007; Richardson, Flollinger, et al. 
2009; Sonnentag et al. 2012; Hufkens et al. 2012). The main species in the vicinity of the 
tower are red maple, sugar maple, and American beech. Qualitative information about the 
autumn foliage colour change stage is summarized in Table 1, and is based on descriptions 
from ~600 local spotters for the foliage network (http://www.foliagenetwork.net). 

Each of these three sites were selected because of the availability of key field data 
needed to assess the utility of the DB phenology product, as well as the sites’ specific and 
fairly strong phenology signal. Following the objective of the article, we use these sites to 
develop a ‘proof of concept’ for using the DB data for estimating phenology parameters 
(noting that further work is needed to test and extend the approach to other areas and over 
time). 


2.2. Daily MODIS VI 

Both the global MODIS operational and daily DB algorithms utilize the semi-empirical 
kernel-driven RossThickLiSparse-Reciprocal (RTLSR) model (Roujean, Leroy, and 
Deschamps 1992; Lucht, Schaaf, and Strahler 2000) to capture the radiative transfer (Ross 
1981) scattering and shadowing components (Li and Strahler 1992) of the surface to recon- 
struct the BRDF at each pixel. This study utilized a daily rolling strategy over a 16 day 
period, with an emphasis on the most recent high-quality observation in the dynamic input 
stream by dropping the oldest and adding the newest looks (Lucht and Lewis 2000; Roy 
et al. 2006; Shuai et al. 2008; Shuai 2010; Ju et al. 2010) to retrieve BRDFs and spectral 
albedo quantities from high-quality, multi-spectral, cloud-free, atmospherically corrected 
surface reflectances. Less change will be found on the retrieved BRDF shape if insufficient 
clear observations are added in the dynamic input stream, especially during the coming 


Table 1. The field-based phenological stages for almond orchard (upper left), forest (lower left), and winter wheat (right). 
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continuous cloud days. A full BRDF retrieval was attempted if at least seven cloud-free 
observations were available during a 16 day period, and these observations adequately sam- 
pled the viewing geometry (Lucht, Schaaf, and Strahler 2000; Shuai et al. 2008). Otherwise 
a poorer-quality magnitude inversion was performed using a priori BRDF information. 
Validation indicates that both full and backup magnitude inversion have a good agreement 
with the available ground sites (Shuai 2010; Wang et al. 2012). To remove the view angle 
effect, the NBAR data were normalized into nadir view with the sun at local solar noon in 
terms of the most recent day of the retrieval period. 

We generated a time series of daily NBAR at the local solar noon for each day in the 
growth period of three vegetation sites. Satellite-based Vis, which provide an indication of 
canopy greenness, are commonly used to monitor seasonal, inter-annual, and long-term 
variations in vegetation structural, phenological, and biophysical parameters (Morisette 
et al. 2009). Vis, spectral transformations of two or more bands, enhance the contribution 
of vegetation properties and allow reliable spatial and temporal inter-comparisons of vege- 
tation activities, and the normalized difference vegetation index (NDVI, Equation (1)) and 
enhanced vegetation index (EVI, Equation (2)) are most widely used in vegetation studies 
(Reed et al. 1994; Myneni et al. 1995; Huete et al. 2002; White et al. 2009; Zhang et al. 
2003; Zhang and Goldberg 2010). As values of NDVI and EVI normally vary with sun and 
view geometry (Moody and Strahler 1994) in the order of 0.05-0.2 (Leroy and Hautecoeur 
1999), it is desirable to use view angle-corrected reflectances such as NBAR for consistent 
comparisons. Therefore, time series of NDVI and EVI for each study site were calculated 
from the resultant daily MODIS 500 m NBAR for each day in the growth period of winter 
wheat (2005-060 to 2005-190), the developing bloom period of the almond orchard (2006- 
001 to 2006-160), and the senescence periods of 2005 and 2006 over the Bartlett forest 
region (from days 195 to 299). 


EVI = G 


NDVI= PNIR ~ Pred , 

PNIR + Pred 

(1) 

„ PNIR ~ Pred 

(2) 

PNIR + ClPred — CSpblue + T 


In the above, p represents the nadir view angle-corrected surface reflectance in the MODIS 
near-infrared (NIR, 841-876 run), red (620-670 nm), and blue (459^479 nm) bands, L is 
the adjustment for canopy background, and C\ and C 2 are the coefficients of the aerosol 
resistance term. The coefficients adopted in the EVI algorithm are L = 1, C\ = 6, C 2 = 
7.5, and G (gain factor) = 2.5 (Huete et al. 2002). 


3. Results and discussion 

3.1. Bloom detection for the almond orchard 

The DB algorithm was used to compute the daily MODIS NBAR-derived Vis (EVI and 
VI) from the 500 m spectral NBAR over the California almond orchard. The NBAR- VI and 
NBARs (NIR, red, and green bands) were extracted for one 500 m MODIS pixel (Figure 1 ), 
centred over a 30 km by 30 km region (NBAR-EVI, Figure 2), to capture the development 
of the almond growth stages (Figure 3) over days 1-160 of the 2006 growing season. Both 
NBAR-NDVI and NBAR-EVI in Figure 1 capture the slow increase in the early growing 
season (days 20 to ~40) as the green tips increase, with a significant decrease in NDVI (and 
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Figure 1. Time series of the 500 m MODIS NBAR-based EVI (cyan) and NDVI (grey) extracted 
from full inversion (solid dot) and backup magnitude inversions (empty dot), the standard 8 day 
MODIS NBAR-derived NDVI (green diamond) and EVI (brown triangle), and spectral nadir BRDF 
adjusted reflectances (NBAR) in red, blue, green, and NIR bands over the almond orchard in the 
northern San Joaquin Valley region of central California with embedded phenological phases (shaded 
with watermark backgrounds) and accompanying pictures. 


a shallow dip in EVI) occurring during the popcorn bloom-petal fall stage (days ~40-70), 
followed by a rapid and continuing increase in the jacket stage and the early nut devel- 
opment stage (days 70-110) and a shift to saturated and stable values through the late 
nut development stage (days ~1 10-160). The NDVI varies from 0.42 to ~0.68, while the 
EVI plot, constrained by the blue band, has a lower dynamic range from 0.2 to ~0.45. 
An increase in understorey in the early growth period, especially in the warm San Joaquin 
Valley, may contribute more to the increase in NIR and VI values than the development of 
green tips from the almond trees. A suppression of NIR and a slight increase in red and 
blue bands coincide with the emergence of pinkish white blossoms (Figure 3). Thus both 
the NDVI and EVI plots record dips as the blooming stage dominates the field of view of 
the MODIS sensor, and dip further when the green understorey is obscured by the fallen 
petals. 

The daily MODIS NBAR-EVI regional maps (Figure 2), also calculated by the DB 
algorithm for days 1-160 of 2006, capture the phenological transition of the larger 30 km 
by 30 km area centred on the same almond orchard in the northern San Joaquin Valley 
region of central California. The transition of EVI values in the centre of the region repli- 
cates the above pixel analysis, with values varying from ~0.2 to 0.4 from dormant to nut 
development stage. Two irrigation ditches, winding through the centre region, may how- 
ever reduce reflectance values of the underlying 500 m pixels. The lower EVI values in the 
upper left of the region are related to a mountainous area (to the west of route 5-West Side 
Fwy) that is covered by very sparse herbaceous vegetation. Relatively high EVI values (up 
to 0.8) are acquired in the lower right of the region, which is covered by a mature almond 
orchard with the rich understorey of pasture, hay, and cultivated crops in terms of the 30 m 
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Figure 2. Temporal EVI maps of a 30 km by 30 km region centred on the investigated almond 
orchard in the northern San Joaquin Valley region of central California for the year 2006, with a 
mountainous region covered by sparse vegetation at the upper left, various densities of orchards in 
the middle, and pasture/hay/cultivated crops (in terms of the 30 m NLCD 2006 classification map) 
at the lower right. 

national land cover data (NLCD classification map). These EVI values drop quickly once 
the understorey green vegetation developed into the senescent period or got harvested after 
day 130. 


3.2. Investigation of the growth phases for winter wheat crops 

Figure 4 shows the seasonal winter wheat phenology profiles from the daily 500 m DB 
MODIS NBAR-derived NDVI and EVI. The dynamic range of the Yucheng site varies 
from 0.25 to nearly 0.60 for NDVI, and from 0.15 to ~0.45 for EVI. Comparisons with 
the field-based record of phenology events (Table 1) show that both the MODIS daily DB 
NBAR-NDVI and EVI captured the main growth features of the winter wheat crop with 
a rapid green-up from the middle of March to early April, followed by a slower increase 
after the emergence of the first node. Then the indices reach their highest values during the 
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Green tip 


Green tip and pink bud stage 
Full bloom 


Popcorn 


Popcorn, bloom, and petal fall stage 


Jacket 


Out of jacket 


Jacket and 
early nuts stage 


Jacket and early nuts developing stage 


Nuts developing stage 


Figure 3. The growth stages of California almond trees from the dormant stage (upper left) to nut 
development stage (lower right); pictures provided by UC Statewide IPM Project, 2000 Regents, 
University of California, http://www.ipm.ucdavis.edu/PMG/C003/m003bcwhybloom.html, and http: 
//bluediamond.com/applications/in-the-field. 


2 week period of heading and flowering, with a gradual decrease during the ripening period, 
and a final drop to the harvest season in late June. The daily NBAR Vis in Figure 4 capture 
some fine details during the periods of rapid change in this growing season. For instance, 
in the flowering and ripening periods, which were captured with a large number of high- 
quality full retrievals, the daily VI plots exhibit a shallow dip, particularly in the NBAR-EVI 
plot. In general wheat flowers are light yellow (see the embedded picture E in Figure 4), 
which may partly decrease the ‘greenness’ of the crop, reducing the reflectance in the NIR 
band and inducing a slight decrease in EVI. When the flowering period is completely over, 
the NBAR-EVI values rebound slightly. 

Figure 5 shows the spatio-temporal variability of daily MOD1S NBAR-EVI from days 
2005-060 to 2005-191 for a 30 km by 30 km area centred on the Yucheng experimental 
site. Generally, the time series of the EVI pictures replicate the winter wheat pixel analysis 
provided in Figure 4, with values varying from ~0.2 to 0.6. The pure white areas in days 
60-62 are affected by cloud, and the cluster of pixels with an EVI <0.2 located in the upper 
left is related to the town of Yucheng. These regional images indicate the applicability of 
field phenology information to the larger local area. Widespread greening occurs across the 
region between days 79 and 87, while widespread maturity appears from day 120 onwards. 
Increasing regional senescence was evident from day 143 onwards, and harvest appears to 
have taken place from day 165 onwards. 
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Figure 4. Time series of the 500 m daily MODIS NBAR-derived EVI (cyan) and NDVI (purple) 
extracted from full inversion (solid dot) and backup magnitude inversions (empty dot), the standard 
8 day MODIS NBAR-derived NDVI (green diamond) and EVI (brown triangle), and spectral nadir 
BRDF adjusted reflectance (NBAR) in red, blue, green, and NIR bands over winter wheat crops at the 
Yucheng site with embedded phenological phases (shaded with watermark backgrounds) and accom- 
panying pictures ( a ) turning green, ( b ) erect growth, (c) first node visible, (d) boot stage, (e) flowering 
and ripening, (f) ripe period, (gj kernel hard period, and (h) harvest. 


* 00*040 

* 00 * 0*1 

* 00 * 0 ** 

* 00 * 0*1 

* 00 * 0*4 

* 00*041 

* 00*044 

* 00*047 

* 00 * 04 * 

* 001 - 04 * 

* 00*070 

* 00*071 

* 00 * 07 * 

* 00*071 

* 00*074 
















* 00*071 

* 00 * 07 * 

* 00*077 

* 00 * 07 * 

* 00 * 07 * 

* 00 * 0*0 

* 0010*1 

* 00104 * 

* 00 * 0*1 

* 001 - 0*4 

*001041 

* 00 * 0*4 

* 00 * 0*7 

* 00 * 0 ** 

* 00 * 0*4 















*00 V 0*0 

* 00 * 0*1 

* 00 * 0 ** 

* 00 * 0*1 

* 00 * 0*4 

* 00 * 0*1 

* 00 * 0*4 

* 00 * 0*7 

* 00 * 0 ** 

* 00 * 0*0 

* 00*100 

* 00*101 

* 001 - 10 * 

* 00*101 

* 00*104 

*00*101 

* 00 * 10 * 

* 00*107 

* 00 * 10 * 

1 

* 001-110 

*001 111 

* 00 * 11 * 

* 00*111 

* 00*114 

i • 111 

* 00*114 

*001 117 

* 00 * 11 * 

*001 114 

* 00 * 1*0 

ijn 

* 00 * 1*1 

HIH 

* 00 * 1 ** 

mmm 

* 00 * 1*1 

kJHB 

* 00 * 1*4 

* 00 * 1*1 

* 00 * 1*4 

VJ 1 

-- _ . 

* 00 * 1*7 

* 00 * 1 ** 

* 00 * 1 ** 

* 00*110 

* 00*111 

* 00 * 11 * 

* 00*111 

* 00*114 

*001 111 

* 00 * 11 * 

ya| 

* 00*117 * 00 * 111 

* 00 * 11 * 

*001 140 

* 00*141 

*001 14 * 

UBI 

* 00*141 

*001 1*4 

*001 141 

*001 144 

* 00*147 

*001 146 

*001 144 

* 00 * 1*0 

* 00 * 1*1 

* 0 O* l*i 

* 00 * 1*1 

* 00 * 114 

*001 111 

*001 114 

* 001- 117 

* 0 W^ 

* 00 * 11 * 

* 001-140 

*001 141 

*001 14 * 

* 001-141 

*001 144 

* 00 *- 10 * 

* 00 * 1 «* 

* 00 * 1*7 

* 00 *!«* 

* 00 * 14 * 

* 00*170 

* 00*171 

* 00 * 17 * 

* 00*171 

* 00*174 

* 00*171 

* 00*174 

* 00*177 

* 00 * 17 * 

* 00*174 
















*oovuo 

* 00 * 1*1 

* 0 O*l»* 

* 00 * 1*1 

* 00 * 1*4 

* 00 * 1*1 

* 00 * 1*4 

* 00 * 1*7 

* 00 * 1 »* 

* 00 * 1 ** 

* 00 * 1*0 

* 00 * 1*1 




















O O.l 0.2 0.3 0.4 0.5 0.4 


Figure 5. Temporal EVI maps of a 30 km by 30 km region centred on the winter wheat field at 
YCES during the above-ground growth period for the year 2005. 
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3.3. Changes in forest senescence 

The hardwood forests of northern New Hampshire (and particularly the area near Bartlett, 
NH) are justly famous for their annual display of autumn colours, with bright yellows and 
vibrant reds. Figure 6 displays the temporal changes of MODIS NBAR-derived EVI and 
NDVI over the Bartlett experimental site, and spans the period from vegetation maturity 
(days 200-220, picture A) through senescence. Senescence, as shown by the EVI plots, 
started in late August, with a gradually continuous drop, while the diminution of NDVI 
started in early September. This may be caused by EVI’s increased sensitivity to foliage 
information in dense forest canopies through a de-coupling of the canopy background sig- 
nal and a reduction in atmospheric and soil reflectance influence (Huete et al. 2002), as the 
leaves of some forest species are turning more colourful by the end of August. The entire 
region near the vicinity of the Bartlett tower has undergone senescence by late October 
(day ~304) as shown in the landscape phenology pictures D and E. Compared with 2005, 
the autumn senescence in 2006 appears to have begun almost 10 days earlier (day ~264). 
In 2006, a very strong cold front moved through the area around day 265, bringing heavy 
rain and up to 31 ms -1 strong winds, which broke some branches and trees, and hastened 
the autumn senescence process of this deciduous forest. Note that the lowest VI values are 
still higher than 0.5 for NDVI and 0.2 for EVI, as a parcel of predominantly evergreen trees 
to the far north of the tower (shown in picture B-E) does fall within the MODIS 500 m 
pixel bin and is also contributing to the VI signal. 
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Figure 6. Time series of the 500 m MODIS NBAR-derived NDVI (grey for 2005 and blue for 2006), 
EVI (red for 2005 and green for 2006) extracted from full (solid dot) and backup magnitude inversions 
(empty dot), and the standard 8 day NBAR-derived vegetation indices (stars) over the mixed forest at 
Bartlett site with the embedded ground landscape phenology photos taken by Bartlett webcam (lower 
pictures A-E). 
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Figure 7 shows the temporal 30 km by 30 km regional subsets of EVI extracted from 
MODIS-derived NBAR at local solar noon and centred on the Bartlett tower location from 
days 195 to 299 for the years 2005 (upper picture) and 2006 (lower picture). The daily 
regional NBAR-EVIs are used to track the transition of the photosynthetically active green 
deciduous foliage from the thick canopy cover (high EVI values) into the vivid reds, yel- 
lows, and browns of autumn senescence. Again it is clear that the senescence in the region 
occurred quite a few days later in 2005 than in 2006. This is further verified by the report 
from the foliage network (http://www.foliagenetwork.net), which indicated that the foliage 
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Figure 7. Temporal EVI of the small area (30 km by 30 km) centred on the Bartlett site for the year 
2005 (pictures above the dark line) and for the year 2006 (pictures under the dark line). 
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Figure 8. Ground foliage colour archived reports over the Northeast US on 2 1 September 2005 (left) 
and 20 September 2006 (right), reproduced with permission from http://www.foliagenetwork.net. 


season of year 2005 progressed more slowly than year 2006, especially in northern New 
England (Figure 8). Warmer temperatures and a lack of summer rain in 2005 delayed the 
onset of fall colour (http://www.wunderground.com) to about day 264, with low colour 
beginning around day 274, peak foliage occurring around day 288, and the region laps- 
ing into dormancy or past peak by day 299. These transition dates occurred about 10 days 
earlier during the 2006 foliage season. 


3.4. Capturing subtle details by monitoring the daily NBAR-EVI 
The detailed signal of the daily NBAR-based EV1 plots reveals the potential of MODIS 
500 m DB products to capture the natural transition of vegetation phenological events, 
especially for plants that have short phenological cycles, such as crops and some orchards. 
Clearly, subtle increases and decreases are revealed in the daily plots for the winter wheat 
and the almond orchards. For instance, light pink almond flowers quickly developed in the 
‘popcorn-bloom and petal fall’ stage and temporarily shaded and retarded the signal of the 
green leaves. Elowever, this induced only a 2.5% relative decrease in the EVI value during a 
1 month period. Therefore, the absolute EVI change (~0.005) of this period is fairly stable 
as compared with the more impressive total absolute EVI change (0.224) from dormant 
to late nut development stage (Table 2). Such subtle variations are only revealed when the 
daily temporal signature is monitored. 

Table 2 shows the length and maximum and minimum EVI values of each individual 
growth phase (determined by the fieldwork) for the three case study sites, as well as the 
relative variation in EVI expressed as a percentage. Various rates of EVI change are driven 
by different plant varieties, by different growth periods, or by different growth cycles for 
a given plant variety (Table 2). Compared with the tiny 2.5% relative change over 25 days 
during the ‘popcorn bloom-petal fall’ stage for the almond orchard, the jacket stage gained 
a dramatic 57.9% relative EVI change in a 36 day period along with the rapid development 
of green leaves and nuts (Table 2(A)). While the winter wheat at the Yucheng site went 
through its entire growth cycle from emergence to harvest within only about 130 days, 
with the fastest EVI increase occurring during the approximately 12 day ‘1st node visible’ 
period and the greatest decrease occurring during the 2 week ‘milky ripe’ period, accord- 
ing to the foliage reports for the Northeast US in years 2005 and 2006, the major EVI 
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Table 2. Statistics of each growth period from the changes in the NBAR-based EVI for the 
California orchard (A), Yucheng winter wheat (B), and Bartlett mixed forest (C and D). 


Study site 

Growth phase 

Days 

Maximum 

value 

Minimum 

value 

Relative 

variation* 

Variation of EVI for 

Dormant stage 

29 

0.2347 

0.2080 

11.8% 

California almond 
orchard during 

Green tip and pink bud 
stage 

13 

0.2449 

0.2131 

14.1% 

2006-001-160 (A) 

Popcorn, bloom, and petal 
fall 

25 

0.2467 

0.2411 

2.5% 


Jacket stage 

36 

0.3828 

0.2523 

57.9% 


Late nuts developing stage 

57 

0.4336 

0.3682 

29.0% 

Variation of EVI values for 

Tiller 

15 

0.1536 

0.1434 

3.1% 

Yucheng winter wheat 

Turning green 

7 

0.1609 

0.1535 

2.3% 

during 

Erect growth 

16 

0.2508 

0.1611 

27.4% 

2006-060-180 (B) 

First node visible 

12 

0.3872 

0.2458 

43.2% 


Boot stage 

13 

0.4229 

0.3787 

13.5% 


Heading 

5 

0.4618 

0.4229 

11.9% 


Beginning flower 

5 

0.4640 

0.4554 

2.6% 


Ripening 

10 

0.4705 

0.4496 

6.4% 


Milky ripe 

15 

0.4406 

0.3181 

37.5% 


Mealy ripe 

6 

0.3182 

0.2612 

17.4% 


Kernel hard 

3 

0.2565 

0.2374 

5.8% 


Harvest 

14 

0.2374 

0.1905 

14.3% 

Variation of EVI for 
Bartlett mixed forest 

Maturity 

65 

0.6585 

0.5344 

24.4% 

during 

Senescence 

26 

0.5326 

0.3294 

63.7% 

2005-200-320 (C) 

Dormancy 

30 

0.3284 

0.1496 

9.3% 

Variation of EVI for 

Maturity 

38 

0.6650 

0.5404 

27.2% 

Bartlett mixed forest 

Senescence 

58 

0.5516 

0.2275 

70.8% 

during 

2006-200-320 (D) 

Dormancy 

25 

0.2542 

0.2071 

10.3% 


Note: * Relative variation means the relative percent of EVI changes for each individual growth stage to the 
change of EVI values in the entire study period for the target plant. 


decreases occurred during the senescence stage, with 63.7% and 70.8% relative change, 
respectively, in the foliage colour change period from days 200 to 320. But, as revealed by 
the daily NBAR-EVI, the start of senescence in year 2006 occurred much earlier than in 
2005 because of the weather, temperature, and rain stress conditions. 

The Collection V005 operational NBAR product (MCD43A4) is only retrieved once 
every 8 days due to archival constraints. While this frequency can usually be used to suc- 
cessfully capture general phenological events, it lacks the temporal frequency to monitor 
these more subtle variations. Furthermore, daily retrievals (which still use multi-temporal 
observations to inform the retrieval but emphasize the observations from a single date), 
may provide additional monitoring information by occasionally obtaining glimpses of the 
vegetation growth even during cloudier periods of the year. When a piecewise logistic func- 
tion of time (Zhang et al. 2003) is fitted via non-linear least squares to the 8 day operational 
MODIS NBAR-EVI data available for each of the case studies, the first-order trend of 
EVI change and the major transition dates are generally captured. However, various subtle 
details in some special growth stages may be missed or smoothed out in the curve-fitting 
process (Figure 9). For instance, intrinsic photosynthetic activities such as the blooming 
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• Daily retrieval — — Logistic fitting O Transition dates detected by logistic fitting 

Figure 9. Comparison of daily NBAR-based EVI (grey dots) with the daily-EVI logistic fitted curve 
(dark-dashed lines) and the detected phenophase transition dates (empty circle) for 2006 California 
almond orchard (a), 2005 Yucheng winter wheat (b), and 2005 and 2006 Bartlett mixed forest 
(c) and (d). 


period (days 40-60) of the almond orchard, the flowering period of the winter wheat, and 
the multiple local peaks in maturity (possibly induced by the variety of vegetation species) 
in the mixed forest are not captured by the logistic trends. Furthermore, the difference in 
timing of senescence between 2005 and 2006 is somewhat masked by the logistic fits. 


4. Summary and conclusion 

This study investigates the ability to detect the natural variation of vegetation phenologi- 
cal events, using the daily MODIS 500 m reflectance anisotropy products generated by the 
daily MODIS DB algorithm. The daily MODIS 500 m reflectance anisotropy products are 
employed to retrieve daily Vis for winter wheat in northeast China, for an almond orchard in 
California, and for a deciduous forest region in Bartlett, New Hampshire. Field-based mea- 
surements, including frequent field specification of surface phenology information at the 
Yucheng experimental site and the California almond orchard, as well as subjective esti- 
mates from local citizen scientists and webcam photos of Bartlett Forest, were acquired. 
Compared with these ground records, the VI trajectories derived from the daily MODIS 
NBAR products capture not only the detailed footprint and principal attributes of the 
phenological events (such as flowering, blooming) but also the substantial inter-annual vari- 
ability. These studies indicate the potential for real-time monitoring of phenological events. 
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especially at the regional spatial scale. The MODIS instruments onboard the Terra and 
Aqua platforms provide a DB downlink for real-time and rapid-response monitoring. This 
daily NBAR product was developed for the DB community and, along with other MODIS 
LIB and L2G science production software, is distributed via the IMAPP (International 
MODIS/AIRS Processing Package) available from http://cimss.ssec.wisc.edu/imapp/. 

Through these daily DB products, the link between traditional ground-based measure- 
ments and satellite-based retrievals can be further investigated to improve our understand- 
ing of phenological phase and the biological life-cycle events in vegetation. Note that any 
satellite-based phenology algorithm performs over the pixel scale with a mixture of various 
biomes, while the ground measurements for particular individual vegetation species often 
record specific phenology events, such as the emergence of crops, bud burst, or flower- 
ing. Thus, future research must be focused on understanding exactly how landscape-scale 
phenological events dominate the signal in moderate-resolution satellite pixels (Liang and 
Schwartz 2009). Furthermore, the possible presence of residual aerosol, cloud contamina- 
tion, or snow will induce additional variation to the satellite-based observations and must 
be acknowledged. However, despite these caveats, these case studies demonstrate the util- 
ity of the daily MODIS DB albedo, NBAR, and BRDF products. Therefore, given these 
results, the 500 m MODIS Collection V006 reprocessing will also be accomplished on a 
daily basis, improving the global community’s efforts to monitor vegetation phenology. The 
detailed work here focused on specific locations and specific dates mainly to demonstrate 
the concept and address the main objective of the article to examine the utility of the daily 
MODIS 500 m reflectance anisotropy DB product for monitoring the evolution of vegeta- 
tion phenological trends. Further work could include extending the use of the DB product 
for regional and global analysis for the entire time series of MODIS data. 
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